This Page Is Inserted by IFW Operations 
and is not a part of the Official Record 



BEST AVAILABLE IMAGES 

Defective images within this document are accurate representations of 
the original documents submitted by the applicant. 

Defects in the images may include (but are not hmited to): 

• BLACK BORDERS 

• TEXT CUT OFP AT TOP. BOTTOM OR SIDES 

• FADED TEXT 

• ILLEGIBLE TEXT 

• SK.EWED/SLANTED IMAGES 

• COLORED PMOTOS 

• BLACK OR VERY BLACK AND WHITE DARK PHOTOS 

• GRAY SCALE DOCUMENTS 



IMAGES ARE BEST AVAILABLE COPY. 



As rescanning documents will not correct images, 
please do not report the images to the 
Image Problem Mailbox. 



I 



4 



THIS PAGE BLANK (USPTO) 



PCI/DK 00/001 40 



S>VK c>o)>(L,C> 




REC'D 


1 S APR 2000 


WlPO 





Kongeriget Danmark 



Patent application No.: PA 1 999 00749 

Date of filing: 27 May 1 999 

Applicant: Leif 0stergaard Harees 

Jacob Knudsens Vej 5 
DK-8230 Aabyh0j 

This is to certify the correctness of the following information: 

The attached photocopy is a true copy of the following document: 

The specification, claims and drawings as filed with the application 
on the filing date indicated above. 





Patent- og 
Varemaerkestyrelsen 

Erhvervsministeriet 



Karin Schlichting 
Head Clerk 




27/05 '99 15:18 FAZ 



61003 



Plougmann, Vinc 



TOFT & PARTN ERS 



Modtaget PD 
2 7 MAJ 1999 



27.niai 1999 

Ny dansk patentansegning 
Lcif Ostergaard Harecs 

Method for detennining bacmodynamic indices by use of toraographic data 
Vor rcf: 22872 DKl 



Plougmann, Vlngtoft & Partners 



C OPENHACtN OfnOL 

POST OFHCEBOX 3007 
DkTioJI COWNHACtN K 



A/$ ASC. NO. 22379S 



m3»H0NI -rfS 35 63 93 00 
TEUfAX 33 63 96 DO 




AAHHUS OFHCt 

POST OFn CE BOX ^9 
DK . H100 AAI WUS C 
mEPHONE -45 »7 M KOO 

TtlEFAX ^ llii »6 DO 

email pv(!>pv^k 

Olt Pknigrrwin 
KnudEffliVInntoll 
Anpr Pon 

Hennk lUstnip Andcrstn 

leirs^ 
Anne Schoubot 
Htniy Scgwrtf 
Marianne JohnnStA 

Mich»4d cwmum 
Ccrl H«y Ittobstn 
AnrwTWttt BlC»"wnn 
AnnoMuie UdMUrwt 

5Lisi« P. Arntien 
\u\ SimorutP 



IngcUborfiw 
K>nru Wigu 
}A6b rid* f^rtdtHben 

Ki-n Vtfi^ncr 

>nt j0^tn Schmidt 
3U*n Maibon 

Mort'n H(UK»^ 

Ar^a Cr&ibnim 

H<idt mcntn 

Nicolas Kiogh 
HnniJk Bagger Olian 
Peror LJudnip 



Tonrt rnU-Mndeelscn 
SupiMn K Addnwni Boston 



Manager ol ft narttt 
)in Himmer 



Pattnts TrademarkJ De5(gni Copyrtghl Ucemlng Documcnutlon 
European Patent Attom«yi European Trademark Anomeyj 
j Coperthagen Aartiui Munkh ReykJavBt Boston Alicante 

I 

i 



3 





27/05 '99 15:18 FAX 



8)004 



1 



r 



METHOD FOB DETERMINING HAEMODYNAMIC INDICES BV USE OF TOMOGRAPHIC 
DATA 

The present Invention relates to a method for determining haemodynamic indices of an 
5 organ or a part of tissue of a mammal, such as a human being, from a time series of 
tomographic data, especially data obtained by means of Magnetic Resonance Imaging 
(MRI). The method produces maps of distribution of blood transit time as well as other 
indices being significant of the dynamics of the capillary tissue flow from tomographic 
images acquired during rapid bolus injection of a contrast agent or tracer that stays 
10 mainly intravascular In the organ. The haemodynamic Indices may by use of the 

present method be obtained with a spatial resolution equal to the spatiat resolution of 
the tomographic data. 

The method comprises the steps of conversion of tomographic Images into data 
15 representing concentrations of contrast agent as a function of time, and the method 
determines the distribution of flow and mean transit time, cither In absolute values or 
relative to the mean value, whereby the methods facilitates the (t) comparison of 
relative flow distributions to that of other tissue (in particular cerebral tissue), such as 
normal tissue, or (ii) quantification of absolute flow distributions In terms of the 
20 associated extraction of a substance with known capillary permeability. 

The use of the methods for examination and monitoring of cerebral conditions of 
humans are of particular interest, but the method may also be applied to other organS/ 
such as the heart, the liver, kidneys etc., or to other part of tissue or part of organs. 



Microscopic heterogeneity of flow is believed to be a main determinant of how 
efficient the delivery of nutrients and pharmaceuticals to tissue takes place. Especialfy 
30 in diseases where delivery of nutrients such as e.g. oxygen is compromised/ 

determination of flow heterogeneity Is therefore crucial to assess the severity of the 
disease. Such diseases include acute cerebral Ischemia, a frequent cause of death and 
the major cause of adult neurological disability in the Western world. 
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So far, the measurement of flow heterogeneity has been limited to the study to 
superficial vessels In the cortex of anaesthetised anfmals by high-speed intravital 
microscopy. There are no pre-existing tools that allow determination flow 
heterogeneity In deeper structures, or on humans as part of non-invasive, routine 
5 diagnostic procedures. 



The study of the delivery of nutrients to the tissue Is currently done by Positron 
Emission Tomography {PET). Due to the costs and lack of general availabilrty of this 
technique, these studies can, however, not be performed in general patient 
10 management. 

A major limiting factor in the development of new drugs in many diseases Is the cost 
of preclinical and clinical trials to determine beneficial effects of new agents. In acute 
stroke, this is typically done by comparing long-term neurological scores of thousands 

15 of treated and untreated patients. The costs of this work, as well as the total number 
of patients available, therefore limits the rate at which novel drugs become available 
for common use. There exists therefore an urgent need for techniques that in 
individual patients predict the progression of a disease or condition and may be used 
for monitoring of said progression, so that the progression for the individual patient 

20 can be assessed whereby e.g. the efficacy of a drug or a substance can be evaluated 
In details from a much more limited number of patients. Such techniques are not 
currently available. 

With the spatial resolution of some tomographic imaging techniques, the dimensions 
25 of vessels are too small to accurately determine arterial tracer concentration levels 
non-lnvaslvely. Instead, image elements containing partly tissue, partly blood vessels 
must be used In order to characterise the mere shape of the arterial Input curve to the 
tissue. In such cases, absolute values of flow and volume cannot be found, and 
therefore a normalisation routine allowing (i) comparison of serial measurements in a 
30 single subject, comparison among subjects and (iii) absolute quantification the case 
of susceptibility MRl of the brain is necessary as described below. 

The present Invention describes and validates a new, non-invasive method that allows 
assessment of flow heterogeneity on generally available tomographic equipment 
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(Magnetic Resonance (MB). Computed Tomography (CT), PET). Furthermore, the 
technique allows indirect assessment of metabolic parameters. 

Rnally, the technique has high predictive power in terms of disease progression in 
5 cerebral diseases such as Ischemia, thereby providing the means for rapid assessment 
of the efficacy of novel therapies. 

DESCRIPTION OF THE INVENTION 

10 It is an object of the present invention to provide a method for determining 

haemodynamic indices of an organ or part of tissue of a mammal, such as a human 
being from tomographic data. 

By haemodynamic indices is understood such indices as distribution of transit times 
15 and parameter characterising said distribution, quantitative haemodynamic parameters 
obtained, from said distribution such as parameters characterising the deviation of said 
distribution from a reference distribution, as well as other Indices being significant of 
the dynamics of the capillary tissue flow. 

20 It is a further object of the present invention to facilitate the comparison of relative 
flow distributions to a predetermined distribution found for a normal organ (e.g. the 
brain) or for an organ having a recognised disease or condition,- such as having a 
tumour of an Identified tumour type. 

25 It is a still further object of the present invention to provide quantification of absolute 

flow"distributions"ln~terms~of"the"aswciatfed~e>ctr~crctidn~^^^ 

capillary permeability. 

It is a yet still further object of the present invention to provide the above methods so 
^ 30 that they are applicable to human cerebral tissue. 

• Thus, the present invention relates to a method for determining haemodynamic Indices 

of an organ or part of tissue of a mammal Including 
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a) determining a time series of tomographic data pertaining to the organ or part of 
tissue during and after a bolus injection of a tracer dose to said mammal, the tracer 
being substantially intravascular In said tissue^ 

5 b) determining a time series of concentration data being indicative of the 

concentration of the tracer in arteries of the organ or tissue from the time series of 
tomographic data, 

c) determining a residue function of the organ or of the part of tissue by deconvolution 
10 of the time seriee of tomographic data with the time series of concentration data, and 

d) determining a distribution of transit times from the slope of the residue function. 

The residue function may also be understood as a normalised impulse response 
15 function characterising the fraction of the tracer present in the vascular tissue at time 
f after a perfect, infinitely sharp input of a tracer In the feeding vessel. For the 
deconvolution a number of known methods may be applied, such as a Fourier 
transform, Box transform, but preferably the Singular Value Decomposition is used. 

20 The distribution of transit times obtained from this method may be used to determine 
a probability density function <PDF) of a normalised haemodynamic lndex# where the 
index Is normalised by the value of the integral of this index. 

One of the haemodynamic indices obtained from the PDF may be a quantitative 
- 25 haemodynamic parameter, such as a parameter obtained from comparison of the 

determlned PDF-and a prevlously determined reference-PDFreTgT-obtalned^Ty use-of — 

the Kolmogorov Smlrnov test- The comparison between the determined PDF and a 
reference may also include identification of plateau's of the curve of the PDF, 
distribution of area under the curve, area under the curve at one or the other side of a 
30 cut-off value, etc. 

The predetermined reference PDF used for comparison may have been found for a 
normal organ (e.g. the brain) or tissue or for an organ having a recognised disease or 
condition, such as having a tumour of an identified tumour type. 
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Furthermore, the present Invention relates to a method in which the quantitative 
parameter may be obtained by performing the steps of 



5 determining the impulse response function of the organ or of the part of tissue by 
deconvolution of the time series of tomographic data with the time series of 
concentration data, 

determining the relative tissue flow from the Impulse response function of the organ or 
10 of the part of tissue, 

normalrslng said time series of concentration data with the integral of said time series 
of concentration data with resp^ect to time, 

15 determining the normalised relative tissue flow, respectively the normalised blood 
volume, of the organ or part of [tissue by use of the relative tissue flow and the time 
sen'es of normalised concentration data, and 

converting said normalised relative tissue flow« respectively normalised blood volume, 
20 to an absolute value for the tlssUe flow, respectively the blood volume, by means of a 
previously determined conversion factor, 

the quantitative haemodynamid parameter being of metabolic significance and 

r 

determined from the PDF and the absolute tissue flow, respectively the absolute blood 
25 volume. 



The normalised relative tissue flow obtained from this method may be used to 
compare tissue flow determined from successive tomographic scanning of the same 
mammal in order to monitor the development of a condition of an organ or part of 
30 tissue end/or to compare tissue flow of different individuals. The comparison may be 
performed because the data have been normalised which means that they ore readily 
' comparable. ; 

j 

f 
I 

I 
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The blood volume Is obtained by dividing the blood "flow by the mean transit time 
(MTT) or as the area under the tissue concentration curve (or the impulse response 
function) of the first pass of the bolus tracer passage. 

5 It is understood that the tomographic data may be obtained from a variety of 

methods, of which MRI and Computed Tomography (CT) are preferred. However, for 
some applications It may be preferred to use other methods, such as Positron 
Emission Tomography (PET) or Single Photon Emission Computed Tomography 
ISPECT). 

10 

The above method may further comprise a step, wherein the normalised relative tissue 
flow, respectively the blood volume^ is also normalised with the ratio between body 
weight of the individual mammal and the Injected tracer dose. By adding this step, the 
determined relative tissue flow Is furthermore made comparable to determined relative 
15 tissue flow from tomography time series of the same Individual mammal in which the 
amount of the Injected tracer dose have been changed. The relative tissue flow 
determined for one individual mammal may furthermore be compared to relative tissue 
flows determined for other individuals. 

20 The conversion factor may be a factor that is generally applicable for the present 
method to members of a mammalian specie. 

In particular, the conversion factor may be a factor that Is generally applicable for the 
. present method to an organ or tissue of the mammalian specie. 
25 

As a special case, a parameter (E) significant for the local extraction of a substance 

may be determined by a method further comprising the following steps: 

calculating the relative flow heterogeneity {w(f)) as a function of the relative 

flow (f) from the distribution of transit times, 
30 estimating a value (P) for the local capillary permeability, 

estimating a value (S) for the local capillary surface area, 
calculating said parameter (E) as the .Integral value of the relative flow 

heterogeneity (w(f)) multiplied by one minus the natural exponential function of the 

negative ratio between 
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i) the product of the local capillary permeability (P) and the local capillary surface 
area (S), and 

ti) the product of the relative flow (f) and the absolute tissue flow {FJ 
with respect to the relative flow (f). 

5 

The substance in question may, e.g., be oxygen, glucose or another important cellular 
metabolic substance or it may, e.g., be a drug or another substance having a 
therapeutic effect. 

10 Estimated values of the local capillary permeability (P) and the local capillary surface 
area (S) are well-known and may be found from standard works within the subject 
area. 

The tomography data for the method according to the Invention are preferably 
15 obtained by means of magnetic resonance imaging, and the method ts furthermore 
preferably applied to cerebral tissue. For cerebral tissue, it has been found to be 
advantageous to obtain the tomographic data by means of susceptibility contrast 
magnetic resonance imaging. 

20 The predetermined conversion factor may be a constant factor applicable for the 
present method for any organ or any part of tissue of the mammalian specie. 
Preferably, the predetermined conversion factor is a constant factor applicable for all 
of cerebral tissue of the mammalian specie, which is strongly Indicated from empirical 
evidence. 

25 

The-tomographlc-data-discussed-8bove-and"used"In~the present method will normally 

comprise information pertaining to subreglons of sections of the organ or part of 
tissue and the haemodynamic indices are determined for at least a substantial part of 
said subreglons. The data of the subreglons may be represented graphically as pixel 

30 values (as a gray-scale image or colour Image of the section), and the quantitative 
haemodynamic parameter determined from the tomograpic data may also be 
represented as Images subdivided* into a plurality of pixels each representing a 
quantitative haemodynamic parameter pertaining to one of said subreglons. 

I 

I 

i 
I 

If 



10 



27/05 '9fi 15:21 PAJ 

■ ^ n ®on 



8 

The present invention furthermore relates to a system for processing of time series of 
tomographic data pertaining to an organ or a part of tissue according to one or more 
of the above described methods according to the invention, said system residing on a 
computer having means for producing an output representative of at least some of the 
5 determined haemodynamic Indices. 

The present invention further relates to a method for evaluating the efficacy of a drug 
or a substance on an organ or on a part of tissue of a mammal by means of 
heemodynamic Indices of said organ or of said part of tissue obtained by one or more 
10 of the methods according to the invention as described above. For this use, a system 
residing on a computer as described in the previous paragraph may advantageously be 
applied. 

The present invention further relates to a method for obtaining information of the 
15 likelihood of recovery of an organ or part of tissue in a living mammal upon or during a 
period of insufficient vascular supply of said organ or of said part of tissue in the 
mammal comprising determining haemodynamic Indices by one or more of the 
methods according to the Invention as described above. 

The present invention further relates to a method for obtaining information relevant for 
discrimination between relevant therapy of an organ or part of tissue In a living 
mammal upon or a period of Insufficient vascular supply of said orgar> or of said part 
of tissue in the mammal comprising determining haemodynamic Indices by one or 
more of the methods according to the invention as described above. 

Furthermore, the present Invention relates to the use of information obtained by use of 
one or more of the methods according to the invention as described above for 
preparing a reference table for use in discrimination of a treatment schedule for an 
Individual mammal or group of mammals for which Information have been obtained in 
30 a manner similar to said information. • 

The equipment for performing the tomographic data for use in the method according 
to the present invention Is general standard equipment available to such an extent that 
the method can be applied in routine patient management. In addition, the 



• 20 
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determination may be performed easily and rapidly for healthy as well severely III 
patients such as patients in coma. The method may be used for diagnosis and 
evaluation of possible further progress of disease whereby a suitable treatment 
schedule may be applied. 



5 



The present method is useful In connection with detection of heterogeneity of transit 
times, velocities or flows by means of external detection of tracers in general, 
whether In biological or technical systems. The technique Is applicable to the 
diagnosis and study of a number of diseases* as well as in the development and 
10 subsequent monitoring of therapies In these diseases. Examples Include: 

(I) Diagnosis and treatment of diseases where, due to altered blood supply, tissue flow 
dynamics and thereby flow heterogeneity of an organ Is altered, e.g. myocardial 
infarction, acute stroke, migraine, carotid stenosis, dementia. 



In particular the method may be used in connection with prediction/assessment of 
subsequent tissue damage based on heterogeneity measurements In the acute phase 
of myocardial or cerebral infarction, evaluation of drug efficacy in these diseases 
based on predictions as to tha subsequent progression of the disease by heterogeneity 
20 measurements, and planning treatment based on prior knowledge of flow or transit 
time heterogeneity and the associated risk of amelioration of disease. 

no Diagnosis, study end treatment of diseases that alter or interfere with vascular 
architecture in any organ or pathological tissue, such that normal flow or transit time 
25 heterogeneity is disturbed, e.g.: (a) Angiogencsis (the formation of new vessels) by 

tumors,3vJiere_the_randQmjformatlQn_o,f_lrr.egularjv^es 

heterogeneity of flows due to the passage through irregular vascular paths, (b) 
Vasculopathies, I.e. diseases where normal vascular wall architecture is disrupted, 
altering the passage and thereby flow heterogeneity of blood or tracer molecules, such 
30 as collagen vascular diseases -(Systemic lupus erythematosus, rheumatoid arthritis, 
other connective tissue disorders), vasculitis (giant cell arterirrtis, polymylagia 
rheumatica), micro- and macroangiopdthies in diabetes, angiopathy due to 
hypertension, (c) Changes In vascular architecture due to chronic or degenerative 
processes of organs {fiver, kidneys, heart, brain) in general. 



15 
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In particular the method may be used In connection with assessment of disease stage 
(e.g. tumor grade) based on heteroeeneity measurements, evaluation of drug efficacy 
In these diseases based on quantification of heterogeneity measurements (e.g. anti- 
5 angiogenic treatment), and planning treBtment based on prior knowledge of flow or 
transit time heterogeneity and the associated risk of amelioration of disease. 

In connection with organ transplantations it Is generally of great value to know the 

vascular supply and the capacity of the vascular system of the organ. 

10 ' i 

i 

(iii) The study of extraction of solutes in normal or' diseased tissue as part of diagnosis 
or treatment. The technique can - given the capillary permeability of a particular 
solute - be applied to a given organ, providing the regional extraction of the solute. 
Examples Include the study of oxygen,jglucose or pharmaceutical uptake in the normal 
15 or diseased brain, e.g. acute or chronic ischemia, dementia, Parkinsons disease. 

Accordingly, the Invention relates to a [method for Identifying a condition relating to 

haemodynamic changes In an organ or tissue of a mammal, diagnosis of disease, 

treatment or prevention of disease, monitoring of disease, and for screening of a 

20 pharmaceutlcally active drug. The method involves use of the determination the tissue 

t 

flow as described herein. 

I 

DETAILED DESCRIPTION OF PREFERRED METHODS ACCORDING TO THE 
INVENTION 

25 : !. 

The present method produces niaps 'of normalised relative or absolute tissue blood 
flow, blood volume and absolute; blood rnean transit time from dynamic, tomographic 
images acquired during rapid boFus injection of a contrast agent or tracer that stays 
mainly intravascular In the organ., \ | 

30 i ; ' : 

By 'blood volume' is understood tracer or contrast agent distribution volume per tissue 
volume. If the distribution volum^ is less than that of total blood (e.g. only plasma), rt 
is understood that conversion to absolute blood volume can be achieved by 
knowledge of regional distributioiji volyme characteristics, e.g. haematocrit, or other 
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normalisation by an Independent techniques (e.g. positron emission tomography) as 
presented in this description. Ulcewise,, 'blood flow' refers to flow of tracer or contrast 
agent distribution volume per tissue volume per unit time, with the aforementioned 
method of conversion to blood flow. Fl'oally* 'transit time' refers to the time taken for 
5 tracer or contrast agent to traverse thB Image element, the arrival at which is defined 
from the arterial input measurement. : 

Furthermore, the method determines the d/str/button of flow and mean transit time, 
either In absolute values or relative to the mean value. The methods further allows the 
(i) comparison of relative flow distributlqns to that of normal tissue (here brain) or (ii) 
10 quantification of absolute flow distributions in terms of the associated extraction of a 
substance with known capillary permeability. 

With the resolution of some tpmogrepbic imaging techniques, the dimensions of 
vessels are too small to accurately determine arterial tracer concentration levels non- 
15 Invaslvely. Instead, Image elements cdntaining paaly tissue, partly blood vessels, must 
be used In order to characterise the mere shape of the erterlal Input curve to the 
tissue. In such cases, absolute values^ of flow and volume cannot be found, and 
therefore a normalisation routine allowing (a) comparison of serial measurements In a 
single subject, (b) comparison among subjects and {c) absolute quantification the case 
20 of susceptibility MRI of the brain Is de^grlbed below. 



OVERVIEW i ; 

I- 

25 The techniques consist of the foIlowin]g steps: 

A. Conversion of tomographic Imageaflnto Images representing concentrations of 

contrast agent as a function of time* 

B. Identification of relative or absolute (i.e. in units identical to those of the 
tomographic Images) arterial tracerj concentrations from th© image data. 

30 C. In the case of refattve arterial tracer concentrations, 

normalisation of the arterial Input area to the Injected dose of contrest agent per kg 
body weight. [ ' 

D. Optional correction of tissue traced curves for delays. 

E. Determination of . I * 
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(]} absolute or relative tissue blood flow 

{ID ttssue Impulse response function 
by deconvolutioo of tissue concentration time curves by tlie arterial tracer 
concentration In each image element. 
5 F. Determination of tissue blood vofume by determining the area under the tissue 
first-pass concentration curve. 

G. Determination of tissue meari transit time. 

H. In the case of MR susceptibility contrast imaging using Gd-chelates In brain tissuet 
conversion to absolute blood'flow and blood volume by a pre-determined constant. 

10 I. Determination of the distribution of flow or transit times in each Image clement 

from the residue function (normalised impulse response function) determined in E. 
J. Comparison of distributions of relative flows to a predetermined distribution found 

for a normal organ; here brain. 
K* Quantification of the distribution of flows In terms of the extraction fraction of a 
15 givan solute with specified capillary'permeabillty, in cases where imaging Is 
performed with MR Imaging with microvascular weighting, or microvascular 
volume can otherwise be inferred from total blood volumes. 



1. Conversion of tomographic data into tracer concentration Images 



20 



The method and associated software handles various modalities, depending on 
whether tracer Injection changes signal Intensity from baseline in a linear or 
logarithmic fashion upon tracer arrival. With a specified option, acquired tomographic 
Images during tracer bolus Injection are converted into concentrations as a function of 
25 time, with measurement time points spaced equally in time. 



1.1. In the case of MR images, weighted towards the transverse relaxation times (T2 
or T2*J, Uypically acquired In brain), the tissue concentration as a function of time, 
C(t), is typically obtained by the formula 

30 



Eq. 1 
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where S(to) is the sisnal intensity before 9ontrast injection (formed as the average of 
signal intensities up to tracer arrivaU, S(tjf is the signal Intensity at time t, and TE is 
the echo time used In the sequence. Here, A: is e constant characteristic of the tissue 
and contrast agent. j :^ 

1 .2. In the case of MR images weighted towards the longitudinal relaxation time (T,) 
images (typically acquired for intravascular tracers in the heart), concentrations ©re 
assessed directly by the change In longitudinal relaxation rate, ARi 



10 



Eq. 2 



where is the relaxlvlty of the applied contrast agent. 

'i 

1,3. In the case of Computed Tomography (CT) Images, concentrations are assessed 
' 15 by the change In Image intensity/ measuVed in Hounsfield units, AH 



Eq. 3 



where k is the characteristic X-ray abso/ptlon of the contrast agent. 



20 



2. Identification of Arterial Vessels 



For Identification of arterial vessels in the Image, the algorithm provides the choice of 
first producing an image that guides this process (2. a.), or proceeding directly to 
_25_calculatlon_of_tissueJlQ_w_(2,bJ_; - 



2. a. If so specified by an option; [the coijicentration time curve of each image element 
Is fitted to a gamma varlate function byinonlinear least-squared regression (by means 



of a fletcher-algorithm) over the 



30 (the re-cIrcuIation is specified by'the usj^r, determining 

!! 



lime Injtjerval until visible tracer re-circulation occurs 



(a) Tracer arrival tiij^a. 



(b) Area under first pass (t^is quantity Is proportional to the local blood 



volume). 
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These two quantities are stored as floating point binary Image file, altowing 
visualrsation as an image where each Image pixel corresponds to quantities (a) and (b), 
respectively. An artery Is then visually identified in these images by 
(a) Anatomical location. 
5 lb) Early tracer am'val relative to the arrival In tissue. 

(c) Largie area under first pass {corresponding to a large blood volume and 
thereby a large portion of the vascular volume being contained In a 
voxel). 

10 2,b, By observing images of tracer concentration as a function of time in a tool that 
allows visualisation of the time-course of single pixels, arterial levels can be Identified. 

Upon Identification of the arterial Input function, the corresponding signal intensity 
time curve is fed to the program below in the form of en ASCII file. 

15 

3. Normalisation of Arterial Input to Allow Comparison of Serial Measurements 



In cases where arterial concentrations are not known In the units of tissue 
concentrations, the algorithm includes a method to allow comparison of relative tissue 
20 blood flow and tissue blood volume In serial examinations in the same subject or 
animal. In order to achieve this, the area of the arterial input function (as determined 
by gamma varlate fitting above) is normalised to the Injected dose per body weight 
prior to the deconvolution below. Likewise, subsequent blood volume measurements 
are normalised to the aforementioned, dose-normalised arterial Input area, 

25 

4,~Determination-of-Blood-Volume 



Using known tracer kinetic principles, relative blood volume Is determined as the area 
under the tissue concentration time curve corresponding to the first-pass of the 
30 Injected tracer or contrast agent (See e.g. Figure 2 in Example 1), normalised to the 
arterial Input area determined in 3. The area is determlr^ed both by direct, numerical 
Integration, and by gamma van'ate fitting (See section 2). 
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5. Determination of Blood Flow and Mean Transit Time 
Tissue flow, Fi, is determined as the solution to the equation 
5 C,(0 = F, .C,(0®/?(0 Eq. 4 

where ® denotes convolution, C.lt) ]s the arterial concentration as a function of tima, 
and R(t> is the residue function, describing the fraction of particles still present in the 
vasculature at time t after a unit Impulse input. 
10 Mean Transit time (MTT) is given from the blood volume (V) and F, the central volume 
theorem 

MTT =r-^ Eq. 5 

With tissue and arterial concentrations measured at equidistant time points ta 
15 =t, + At,....tN. Equation 4 can be reformulated as a matrix equation 











0 


0 ^ 
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m 





















£q. 6 



The equation is (after applying Simpsons rule to the matrix elements (0stergaard 
20 1996b}) solved for F|*R(t) in each pixel by Singular Value Decomposition (SVD)^ with a 
cut-off in diagonal eigenvalues at 20% of the maximum value in order to suppress 
experimental noise typical of tomographic Images (signal to noise ratio of the order of 
10). The noise cut-off of 20 % is user specified. In the case of higher slgnal-to-nolse 
ratio ISNR) or specific SNR characteristics, a Monte Carlo simulation scheme is 
25 provided to optimise the cut-off (06tergaard 1996b). 

The algorithm allows Independent correction for tracer arrival delays prior to 

decon volution by shifting the tissue curve in a given pixel by the delay fitted in 2. a** 
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thereby reducing bias fn flow rates due to simple delays of the tissue curve relative to 
the arterial input curve (0stergaard 1996a). 

The algorrthm calculates the followmg quantities, and stores them es digital, floating 
5 point Images: 

(a) Tissue blood volume (As determined section 2. a. and 4). 

(b) Tissue blood volume (As determined by simple, numerical integration 
up to tracer re-circulatlon - section 4). 

(c) Blood mean transit time determined as the area under the residue 
10 function response RCt). 

(d) Relative or absolute tissue flow as determined by the maximum value 
of the impulse response function <F,-R(t|| - Section 5. 

Here, relative flow, F^, means a value proportional to (by the same 
factor throughout the digitised images of relative flows) absolute flow 
15 In cases where absolute arterial flow values cannot be determined. 

Again, absolute flow refers to cases where arterial concentrations are 
known in the units of the tissue concentration data, 

(e) Relative or absolute tissue flow as determined by Equation 5., using 
the blood volume determined in (b> divided by the MTT determined In 

20 (c). Here, relative flow means a value proportional to (by the same 

factor throughout the digitised Images of relative flows) absolute flow 
in cases where absolute arterial flow values cannot be determined. As 
above, an absolute flow refers to the cases where arterial 
concentrations are known in the units of the tissue data. 
25 If) MTT determined as tissue blood volume (b) divided by tissue flow (d) 
(Equation 5>. 

(g) The tissue residue function, R(t), 

(h) The goodness of fit, x*. of Equation 4 to the tissue concentration time 
curve. 

30 

6. Absolute Quantification 
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To allow absolute quantification of cerebral blood flow (CBF; ml/lOOml/min) or 
cerebral blood volume (CBV; ml/joOfnl) measurements with Gd-chelates and dynamic 
T2 weighted susceptibility contra;^ agents values are determined as 
CBF - 0.87- Bo- khuhnans) 



CBV - 0.87- Bo- V,: :[humBns> 



CBF = 1 .09- Bo- F, 



j (porcine 



data) 



CBV - 1 .09- Bo* Vr! '(porcine data) 



Where Bo Is the field strength of thel magnet In Testa, and V, and F, are the relative 
volumes and flows determined asfidescribed In section 4 and 5, respectively, where 
10 Impulse response height and tissJe burve area have been eppropriately normalised to 



the aaerial input area. 



ii 1 



7, Determining the Distribution ol Transit Times 



15 To determine the distribution of Jis^e transit times, h(t), the negative slope of the 
tissue response is determined in |ach point as the average of slopes of straight lines 
connecting the previous end follo^vrng point on the Hit) curve (See Equation 13 and 
Equation 14). i r 



Special cases: 



i I 



20 



7. a. At the Initial pointlof R(t), where only one slope can be defined, the 
slope of the straight Jlirio connecting the initial and second measurement 
of Hit) Is chosen ; r j 

7.b. In cases wherejldue to noise, Bit) between to points deviates from 
being either constani cr decreasing ©s a function of time, the slope Is set 



25 to zero. ! j| 

The resulting curve, with the tlr|\| i 



30 



illustratedln Exampfe~2 (Osterga^rc 

•r 
:( 

8. Determining the Distribution Felatlve Rows. 

II 

In the remainder of the descriptlcyi 



xls normalised to the mean transit time, is 
ig98ar^ 



relative flows refer to a defined fraction of the 
mean flow, i.e. a dimen6ionle88;c[u'antlty, independent of any absolute quantification 
of flow- 
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Under the assumption of equal 
times can be converted into a distribution 
such that the mean of the di 



lengths of capillary paths, the distribution of transit 

w{f) of relative flows f (i.e. flows normalised 
stribufion is 1 ) by the relation 



where V is the associated blood vplume 
to the normalisation In Equaxlon 9 



(Units of V and may here be arbitrary due 
, and thereby 



10 



where t is a given transit time in t 
negative slope of R(t) and dt is 
thereby R{t)). 



given 



15 



Its Is finally required that the 
(PDF), Ke. has area 1, and has 
flow). I.e. 



20 



TTiis IS achieved by Iteration, 
above, in each step scaling f and 



121021 
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Eq. 7 



Eq. 8 



hie distribution above, h{t) the corresponding 

by the time resolution of the measurements (and 



resulting curve w(f) be a probability density function 
mean value 1 (i.e. describes flows relative to total 



to/- Jw(f)df = 1 
0 • 



Eq. 9 



numjerically integrating the expressions of w(f) and f 
lA^(f) by the preceding integral. 



9. Determining Abnormal Flow Distributions 



25 



The resulting distribution of relativ^e 
determine its values at predefinec 
normal tissue (Figure 4b in £xam^li 
The values read 



30 



flows is recalculated by linear interpolation to 
relative flow values previously determined for 
el). • 

TABLE 1 



21 



# 



27/05 '99 15:24 FAX 



In each Image element, the distribution 
5 of a nonparametric Kolmogorov Siinirnov 
calculated. 
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f 

3.8053 


1 

y(0 

L65E-03 


f 

0.6941 


1.4594 


2 a65 


* 0.0226 


0.647 


1.6915 


2.3306 


1 

: 0.0372 


0.6205 


1.6279 


1 .7879 


. 0.0627 
t 


0.5969 


1.5812 


1 .5881 


: 0.1198 


0.5815 


1.351 


1.4459 


' 0.1751 


0.565 


1.3393 


1 .3535 


' 0.2491 


0.5448 


1.2367 


1 .2785 


^ 0.3303 


0.5262 


1.1433 


1.1912 


' 0.4058 


0-5133 


1.0688 


1.1108 


1 0.4793 


0.5007 


0.928 


1.0474 


0.64i8 


0.4822 


0.6559 


0.9972 


• 0.7218 


0.4669 


0.6366 


0.9438 


: 0.6032 


0.4571 


0.4282 


0.8981 


0.986 


0.4416 


0.3133 


0.8214 


i 1.1098 


0.4316 


0.0244 


0.7794 


I 1.3562 


0.42 


0 



Ts compared to this standard curve by means 
(KS) test# and the corresponding p value Is 



The program returns digitized images of 
8. a. The KS probability 
10 normal distribution. 

8^b^20-lmages-cont^ning 

Table 1. 

8.C. With a specified 
of Is performed (with[.the 
15 times or relative tranjfrt 



p as a quantitative measure of deviation from the 
the-values-of-w(f)-at-th8-predefin©d_valuB5JJn- 



'option, ^corresponding comperlson of the distribution 
curye in Figure 4a in Example IJ, and the transit 
i'rt times displayed. 



10. Absolute Quantification of FIciw heterogeneity and Solute Extraction 
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With a specified option, for susceptibility contrast MR imaging, the program combines 
the relative flow heterogeneity w(f) with absolute flow estimates (See Section 5) to 
calculates distribution of absolute flows In ml/1 OOml/mln. Furthermore, using the 
microvascular weighting (For example asidemonstrated in 0stergaard 1998c) and 
5 table values for the microvascular dimensions of the organ in question {for brain 
capillaries e.g. 8 \im diameter and 1 20 ^Tjn length), blood volume is converted into 
capillary surface area, S. The local extraction E of a substance with capillary 
permeability P (as given by literature values). Is then given by: 

-i 

10 E = \wif)^^ -e"^ycff Eq. 1 0 

I 

! 

This quantity Is visualised as a separate image, based on the specified value of P. 
In cases where tomographic images do not possess a specific microvascular 

weighting, S can be inferred from the blood volume by literature values describing 

I 

16 microvascular dimensions and fractions 1|or the relevant organ, 

I 

i 
! 

i 

i 
t 

i 
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EXAMPLE 1 . 

MODELING CEREBRAL BLOOD FLOW AND FLOW HETEROGENErTY FROM MAGNETIC 
RESONANCE RESIDUE DATA 

5 SUMMARY OF THE EXAMPLE | 

Existing model-free approaches to determine cerebral blood flow by external residue 
detection show a marked dependence of flow estimates on tracer arrival delays and 
dispersion. In theory, this dependence can be circumvented by applying a specific 

10 mode? of vascular transport and tissue flow. heterogeneity. A method Is presented to 
determine flow heterogeneity by tflR residue detection of a plasma marker. Probability 
density functi'ons of relative flows measured in 6 healthy volunteers were similar 
among tissue types and volunteers, andjwere In qualitative agreement with literature 
measurements of capillary red blood cell and plasma velocities. Combining the 

15 measured flow distribution with a model of vascular transport yielded excellent model 
fits to experimental residue data. Fitted <gray:whtte flow rate ratios were in good 
agreement with PET literature values, as well as a model-free singular value 
decomposition (SVD) method in the 6anr>e subjects. The vascular model was found 
somewhat sensitive to data noise, but showed far less dependence upon vascular 

20 delay and dispersion than the model-free SVD approach. 

BACKGROUND j 

A technique to determine cerebral blood flow tCDF) by magnetic resonance (MR) bolus 
tracking of an intravascular contrast bolus has recently been presented (0stergaard 
25 1 996a}. By performing non-parametric singular value decomposition (SVD) 

deconvofution of tissue time concentration curves by a non-lnvaslvely determined 



arterial input function^ the algorithm (hereafter referred to as the SVD method) 

generates pixel-by-pixel maps of CBF (fijstergaard 7 996a). The SVD method has been 

demonstrated to yield CBF values in excellent agreement with PET in normal 

J 

30 volunteers (Ostergaard 1998a), ais well as in an animal hypercapnia model (0stergBard 
1998b>. Although the model-free SVD rpethod offers the advantage of being 
independent of the underlying vascular structure, the method Is somewhat susceptible 
to dispersion and delay of the measured AlF be.fore It reaches the imaging pixel 
<0stergaard 1996b). Especially In the setting of major vessel disease* dispersion In the 
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I 

i 



feeding vessel may be significant relative to tissue tracer retention, causing 
underestimation of CBF, and therefore o>!rerestlmation of the CBV:CBF ratio(0stergaard 

1 gS6b). This ratio, the ptasma meen transit time (MTT), Is 9n important parameter in 

i ! 

evaluating cerebrovascular perfusion reserve, and therefore the InabilitY of the SVD 
5 approach to distinguish vascular dSspersion from prolonged tissue MTT may ultimately 
Impair its clinical use. j 

= • ! 

Basslngthwalghte and co-workers! have developed modelling tools to describe major 
vessel transport as well as microvascular tracer retention {King 1993; King 1996). A 
10 derived model of the coronary circulation has been successfully applied to MR data, 
allowing non-Invasive measurements of coronary blood flow (Kroll 1 996). This model 
(hereafter referred to as the vascular mo'del), modified for the cerebral circulation, may 
ultimately provide estimates of c6f andjMTT independent of major vessel delay and 

dispersion. The aim of this study was .ta[ extend this vascular model to the cerebral 

; i 

15 circulation. First, a first-order expression of flow heterogeneity In the cerebrat 
circulation was derived by a model-free ianatysis of tracer retention In areas of 
negligible major vessel dispersion.* To validate the model, the vascular model was 
fitted whh measured flow heterogeneitv^ to human MR residue measurements, and 
compare the resulting flow rates to literature values, as well as the SVD method. 

20 Finally, the sensitivity of the vasiiilar model to tracer delays was compared in 
simulated data. ^ 

i 
I 

THEORY . : 

Vascular model • 
25 A modified version of the vasculajr model previously described by Kroli et a] (Kroll 

1996) was used. The vasculatur^* was modelled a major, feeding artery in series with 
^2G-small vessels In-paralleF- Seel Figure|^ln-the-vascular model the-f low-is directed — 

along the parallel pathways, withjwi representing the fraction of flows with values 

1 • t 

fi'Fp, where Fp is total flow. The feeding artery was described by a fixed, relative 
30 dispersion in transit times (RD = d.4B) apd a delay, determined by Us volume fraction, 
V„ (King 1993). The capillaries vlere modelled as simple delay lines of fixed length 
(100 ^m), with volume Vp, In theifollowjng, the tissue flow Fp, the feeding vascular 
volume V^, and capillary volume jVp was allowed to very. The observed signal 
changes due to magnetic susceptibilrty Contrast agent arises - when using a spin echo 

■■ i 
•• . \ 
1 



25 



27/05 '09 15:26 FAX 



(21026 



; 

sequence (See Materials and Methods below) - mainly from capillaries (Rsb! 1991; 

Boxerman 1995; Weisskoff 1994): In th^ model, total tissue tracer concentrations 

were therefore calculated based on the amount of tracer in the small, parallel vessels. 

The distribution of transit times in:the capillary bed was incorporated by an algorithm 

5 assigning appropriate flows and weights Ito the parallel vascular paths, to achieve a 

given heterogenerty IKIng 1996). This flow heterogeneity is described by a probability 

density function (PDF), assigning d probability w(f) to a given relative flow f , I.e, flow 

relative to the mean flow, Fp. In the follo[wing, it Is described how an estimate of flow 

I 

heterogeneity In humans is obtained frorn MR residue data. For Q more detailed 
10 discussion of modelling flow heterogenit/, see King et ai, (King 1996). 

Row Heterogeneity from Residue Data j 

The fact that the impulse response to a plasma tracer can be estimated by 
nonparametric deconvolutlon of the tissue residue during the tracer passage by a non- 
15 Invasively determined arterial Input furiction (AIR is utilised. From this the distribution 
of plasma transit times, and - under certain assumptions regarding the distribution of 
capillary lengths - the distribution of fltows In the region are derived. 
The tissue concentration Ct(t) of tj-acer In response to an arterial input function C.{t) is 
given in Equation 4. The formula fe equivalent to the Integral 

50 : : 



; • : 0 

» 'i Eq. 11 

where Fx is tissue flow and R Is tfie rBsjifue function. I.e. the fraction of tracer present 
in the vasculature at time t after a perfect, infinitely sharp Input in the feeding vessel. 

Assuming-the-arterial-and-cerebrajTcontpntrations-are-measured-at-equally-cpaced 

25 time-points this equation cah!be discretized, assuming that over short time 

intervals At, the residue function and'artcrial input values are constant in time: 



0 i "}\ '-0 



Eq. 12 
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which Is equal to the expression Equation 6. 

As previously described {.0BterQ^djri 1996b). this equation can be modified to residue 
and arterial input functions that vary tinearlly in tinne. The SVD approach provides a 
5 powerful numerical tool to solve ^quatibn 6 in the presence of experimental noise to 
yield the residue function. The distribution of transit times, h{t), is then found from 

\ .*' 



] Eq. 1 3 

j :i 

10 I.e. the slope of the residue functfon. At a given time point t,. h(t) can be estimated as 

I 

4 

I I- 

The probability density function (PDF);iof transit times (t) Is turned into a dlstrIbutlor» of 
15 relative flow rates f, w{f), by requiring;' 

i ' 

w[f)dfih(t)dt Eq. IB 

] ' 

and that the central volume theorem (Stewart 1894) is obeyed (Equation 7). 
20 Assuming all vascular paths havej equal volume, the distribution of flow rates Is 

obtained from Equation 8, and ndrmairsed to have unit mean flow and area (Equation 
9). I 

i :! 

MATERIALS AND METHODS \ 



25 



Volunteer Data j r 

Six normal volunteers (Age 29 ±4 yrsj were examined according to a standard 
perfusion protocol on a GE Sign»|l,5 ;T imager (General Electric, Waukesha, Wl) 
retrofitted for EPI capabilities {Instascin, Advanced NMR Systems, Wilmington, MA), 
30 using spin echo (SE), echo pianafi imaging (EPI) with a time of repetition (TR) of 1 ,0 



:! 

'I 
I 

il 
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seconds, and a time of echo (TE) of 1C(0 ms. The slice thickness was S mm with an 
In-plane resolution of 1 .6 mm by 1 ,6 rj{m in a 40 by 20 cm fleld-of-view IFOV). A 
total of 52 Images were acquired, starting 15 seconds before l.v. Injection of 0.3 

mmoi/kg contrast medium {DysprosiunS-diamide, Nycomed Imaging, Os!o, Norway). 

'1 

5 Intravascular contrast agent concentrations, C(t), were estimated assuming a linear 
relationship between concentration anjj change In transverse relaxation rate, AR2 
(Villringer 1988; Weisskoff 1994) (Equation 1). 

feeding arterial branches were Identified in the Image slice as pixels displaying early 
concentration increase after contrast liijection (Porkka 1991). This approach does not 
10 determine absolute arterial tracer levels, but provides the shape of the AIF. To 

standardise the analysis below, the arterial Input function was therefore scaled to yield 
a mean CBV of 3% (©stergaard 1998ii). in volunteers, a single arterial input function 
In the imaging plane was used for all flssue regions. 

i 

^1 

15 Determination of Flow Heterogeneity |[ 

Tissue concentration time curves werej formed using Equation 1. Three gray and two 
white maner tissue regions consistingrjof 4 imege pixels (0.05 cc) were then chosen, 
based on cerebral blood volume maps!;(Roscn 1990). The tissue residue function was 
calculated by SVD deconvolution of tijje tissue concentration time curve with the AIF. 

20 The resulting residue function was th^n converted Into a probability density function 
(PDF) of relative flows as described la;the Theory section above 

:i 
.1 

Model Analysis *! 

The experimentally determined flow heterogeneity PDF was then entered into the 
25 vascular model described above. For sixteen gray and white matter regions (0.26-0.4 

cc), Fp, Vp, and V,^ were adjusted to obtain optimal fits to the correspondin g tissue 

_ q 

concentration time curves by nonlfneaj- regression analysis (Chan 1993). The initial 

conditions were Fp«40 ml/1 OOml/mlfij Vp«2%, end V„ = 0.1%. The remaining model 

i! 

parameters are given In Figure 1 . ;l 



30 



i 

Comparison with model-free dpproach| 

To compare the flow rates obtained \|jith the vascular model with those of the model- 
free SVD epproach (0stergaard jet al.l!!l99ea), the height of the deconvolved tissue 
response curve (F,-R(t|) In Equation Stl was determined for the same regions as used 
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I 

I*i 

in the model analysfs above. After dete,rmlning mean white maner flow rate, 9 relative 

|j ' 

gray; white flow ratios were calculated ^or each volunteer, end compared wFth those 
obtained by the vascular model. • \* ' 

I 

5 Sensitivity to Tracer Delays 1.! * 

In volunteer 4, each pixel tissue conce,otratlon time curve was delayed in steps of 
0.25 second by linear Interpolation to ^imufate the effects of tracer arrival delays. The 
simulated Image data were then analysecf as described above, and fitted flow rates by 
the SVD and vascular model approach^plotted as a function of de)ay for comparison. 
10 \' 
Sensitivity to Noise and initial Conditions 

To determine the overall sensitivity of ^parameter estimates with the vascular model to 
experimental noise, a set of synthetic faata sets were generated using the vascular 
model itself, and two sets of typical values for flow, volume, and feeding artery 
15 volume (Fp = 20 ml/100 mf/min, Yp«=2%i and Fp-50 mi/lOO ml/min, 
Vp-3%, V« = 0.5%). These were cortyerted into a MR slgnel intensity time curve 
using Equation 1, to generate a f y pic a^- signal loss (25% for gray matter, equivalent to 
the higher flow) during a bolus passagje. Random gausslan noise was then added, and 
'noisy' concentration time curved werp again calculated from Equation 1 . Simulated 
20 SNRs varied from 400 down to 12, ti|5 latter being typical for raw, pixel-by-plxel data 
obtained with perfusion protocoljs on i clinical MR system. The synthetic curves were 
analysed using the vascular model, u^ng two different sets of initial conditions. These 
were chosen to represent two extrem^'^ of physiological values: Fp-80 ml/100 
ml/mln. Vp-6%, V^ = 1%, and Fp = 20 ml/IOOml/min, Vp-2%, V„ = 0.1%. For each 
25 SNR, 24 simulations were performedjand the mean and standard deviations of the 
fitted model parameters were cafcula&d for further evaluation. The dependency on 



Initial conditions was evaluated by redo^dlng the number of simulations where two 
different initial conditions caused resisting fitted parameters to differ by more than 
1 0% from their mean. To compare stability of the vascular model and SVD method to 
30 experimental noise, flow rates vfere c^termined from the same synthetic curves hy the 
SVD method. 



RESULTS 



I 



( 1 
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Figure 2 shows 9 set of typica! tissue end arterial concentration time curves obtained 
from vorunteer 1. The tissue ROIs consisted of 25-35 pixels (corresponding to 0.3 - 
0.4 cc volumes). The average signal tcj noise ratio (SNR) (defined as the maximum 
tissue R2 increase during bolus passage divided by the standard deviation of the noise 
5 relative to pre-bofus base-line Im^ge iniensity) for the tissue volumes used for model 
validation CO. 25-0.4 cc) was 30.Vhe SN|R of gray matter was a factor of 2 to 3 higher 



than that of white matter, due tof the higher blood volume. 



Flow heterogeneity 



10 Rgure 3 shows the location and size of jhree regions chosen for determination of flow 
heterogeneity in volunteer 4. ThiS regions are overlaid on a CBF map calculated by the 
SVD-method to illustrate the co^rast apd spatial resolution of the techniques. Also 
shown are the corresponding flow heterpgenelty plots derived from the three regions. 
The transit time and derived flo^A> hetejdgeneity PDFs were found to be remarkably 

15 similar among regions and among volupteers. Figure 4a shows all pairs of relative 

transit time t and corresponding^|i(t) m^easured for all regions In all 8 volunteers. Figure 
4b shows - under the assumption of e^iOal capillary lengths - the corresponding plot 
of relative fiow f and w(f) measured for, all regions and volunteers. Due to this 
constancy across regions and subject^ ;the (f,w(f}) points were consequently averaged 

20 into 30 points (Full curve) and u^^ed as |a global expression for flow heterogeneity in 
normal tissue in the subsequentimode analysis. The parameterlsed w(f) is given in 



Table 1. The distribution of fiovsi is rr^rkedly right-skewed, with the majority of 



capillaries having flow rates lessjthan 
reached at roughly 2/3 of the mean fidw. 



25 



he mean flow. The maximum probability is 



Model validation 



In the following, the experimentally d 



ermined flow heterogeneity PDF was applied, 

■I 

described by only three parameters^ V^, Fp and 



and the vascular transport is the|'efore c 

Vp. Figure 5 shows a typical se^of :tis! i^e concentration time curves as well as the fits 
30 provided by the model (Basal gangh'a: 
matter: CBF = 48.6±1.2 ml/1 OOml/rnin 



Notice the model takes into account 
higher flow rates. The quality of mod 



diBF = 60.3±1.0 (SE) ml/IOOml/min, Gray 
white matter: CBF = 19.5±0.8 ml/1 OOml/mln). 
t1 14 observed earlier tracer arrival In tissue with 
eI [Fits to experimental data shown In Figure 5 Is 
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typical for the patients exemined. Table 2 shows the meart grayrwhite matter flow 
ratios for 8 regions. The mean gray:white flow ratio was 2.89±0.35 with the vascular 
modal. 

S TABLE 2, Regional grayrwhite flow ratios in 6 volunteers 

Volunteer no. 



10 



15 



20 



Region 




1 


2 


3 


4 


5 


6 Average 


Basal ganglia 


left 


4.19 


3.41 


3,70 


2.73 


3.76 


3.39 


3.53 




right 


3.76 


3.41 


2.99 


2.80 


3.66 


3.86 


3.41 


Frontal medial 


left 


2.89 


2.91] 


1.94 


2.32 


2.55 


2.03 


2.44 




right 


2.66 


2.73 


3.6 


2.24 


3.06 


3.03 


2.88 


Temporal 


left 


3.55 


3.16 


3.14 


2.55 


4.31 


2.29 


3.16 




right 


4.&2 


3.59 


3.27 


2.02 


3.71 


2.78 


3,31 


Occ.-temporal 


left 


2.50 


2.37 


2.66 


2.20 


2.55 


2.74 


2.50 




right 


2.78 


2.06 


3.70 


2.08 


1.70 


1.97 


2.38 


Occ. medial 


left 


3.60 


2.87 


2.27 


2.06 


2.30 


2.84 


2.66 




right 


2.95 


2.li? 


3.57 


2.17 


2.15 


2.43 


2-57 


Average 




3.34 


2.86 


3.08 


2.31 


2.97 


2.74 


2,89 
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Comparison with model-free approach 



Figure 6 shows gray:whlte matter flow ratios determined by the model approach 
plotted versus corresponding ratios^in identical regions determined by the SVD 
approach. The three volunteers! were chosen based on a significant spread In 
30 individual, regional gray:white mat^r ratios. In order to facilitate comparison between 
approaches, in the figure, linear regression lines for each volunteer are shown 
(Volunteer 1 : y= 1 .02 + 0.21 (r*4o.75). Volunteer 2: y - 1.1 6x - 0.68 (r' -0.91), 
Volunteer 6: y •»0.95 X -0.72 {r* = p.9S)J. The line of identity was within the 95% 
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linjear fix. Notice regional gray: white matter ratios 
the tine of identity- 



and the SVD approach, respectively 
of 90 Image pixels (1.1 cc). Vasculi 
identical regions. Tho SVD approacV 



Sensitivity to tracer arrival delays 

Figure 7a shows the effect of AIF dfilay on fitted flow rates for the vascular model 

The gray and white matter tissue ROIs consisted 
mode! and SVD values were obtained from 
progressively underestimates flow rates with 
tracer arrival delay. The relative undjerestimation is roughly proportional to the delay, 
10 reaching ^5% for white matter andj35% for gray matter at a delay of 3 seconds, 
respectively. For the model approacli, flow estimates are remarkably Independent of 
delay. Figure 7b shows the firted fdpding artery volume as a function of delay. The 
vascular model Interprets increasing delays as increased feeding vessel volume In 
accordance whh the definition of the vascular operator. As vascular dispersion Is a 
15 pWb// unknown in actual measuremlnts, fitting was performed assuming a vascular 
dispersion RD"0.48, although the delay was simulated to be dispcrslon-less Notice 
fluctuations of fitted flow values araund 1 Is accompanied by fluctuations of the fitted 
vascular volumes. These fluctuatiorts and the tendency of flow and vascular volume to 
co-vary are discussed further belov 

20 

Sensitivity to noise and initial cond^ ions 

Figure 8 shows the means and star dard deviations of the fitted flow rates for two 
sets of simulated data (V«-0.5%. Fp=60 mi;i OOml/min, Vp-3%; V«=.0.5%. 
Fp»20 mi/lOOml/min, Vp- 2%), us ng the model (Figure 8a) and SVD (Rgure 8b) 
25 approach, respectively. Raw image data noise was varied from that of typical, clinical 
data 1 2) to 400. The SVD apprc 3ch overestimates flow rates somewhat for this 

choice-of-modeLparamcters^-where fs_the_vascuIar_modeLfitS-are_roughIy_equaLto_tbe 

\ arror bars In Figure 6a and 8b Indicate one 



Input parameters. The uncertainty 
standard deviation, derived from tha simulated data) on fitted flow rates display the 



30 expected increase as a function of 



ncreasing raw image data noise. 



For low SNR, error bars are roughlv 
respectively. For high SNR, however 
reach zero as was the case for the 



equal in size for the SVD and model approach, 
, the uncertainty on vascular model fits did not 
SVD approach. To investigate whether factors 
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Other than nofse contributed to the observed behaviour, the dependence of model fits 
on Initial coriditlons and the tendency of parameters to co-vary was analysed. The 
fraction of fits, where initial conditions were found to significantly affect the fitted 
flow rates {defined as cases where two different initial conditions resulted in fitted 
5 flows that differed by more than 1 0% from their mean) was found to be negligible for 
SNR above 20. For lower SNR, 15-20% of fits gave ambiguous results. Therefore, the 
standard deviations for the low SNR may be somewhat underestimated due to bias by 
the choice of initial conditions. 

During simulations, fitted feeding artery volumes and fitted tissue flows were found to 
10 co-vary: high, fitted flow rates were hence often accompanied by high arterial 

volumes. Figure 9 shows fitted flow rates versus corresponding fined feeding artery 
volumes from simulated curves with SNR = 40. This pattern of co-varying fitted flow 
rates and feeding artery volumes was found at all noise levels. It was found that 
varying flow and feeding artery volume in proportion lead to only small changes in the 
15 shape of the resulting concentration time curve. Therefore, the presence of modest- 
experimental noise leads to relatively large uncertainties In fined flow rates. This is 
thought to explain the unexpected, large standard deviation of flow estimates at high 
SNR for the vascular model, as well as the fluctuations around unit relative flow In 
Figure 7. < 



20 



DISCUSSION 



5 



Overall validity of model 

The vascular model, after incorp({rdtIon of the experimentally determined 
25 heterogeneity PDF, provided excellent fits to the experimental data. The model fits of 
CBF'yleldBd-B-mean~9ray:whlte ffow-ratlo~of~2T8iO:357l^^ 



f few ratlorof ~j 
3cts of similar . 



literature ratio of 2.65 for subjects of similar age (Leenders 1990). The ratios were 

also In good agreement with those found using the model-free approach. In contrast 

to the SVD approach, the model approach provided fits to experimental data that were 

30 essentially Independent of vascular delay. The model was simplified considerably by 

the use of one flow heterogeneity, PDF for all types of tissue. By characterising the 

model by only three parameters, p remarkable stability of the model was obtained, 

even compared to the pixel-by.pi?cel based SVD approach. Co-varying vascular volume 

and flow rate in model fits to lov|r;SNR data was found to be a major contributor to 

>. 
I 



33 



?7/03 98 15:29 FA X ^ ,3,^54 



!31 

t 

uncertainty in model fits of CBF. Below, the IndivicJual elements of the model aro 
discussed in further detail. ^ 

I 

Model considerations i 

5 The model presented here differs slight'ly from that previously described by Kroll eta/. 

(Kroli 1996) in the heart. Kroll et aL (Kroll 1996) described the arteriolar and capillary 

compartments separately^ assigning Tt^Qd relative dispersion and variable volumes to 

the arteriolar vascular paths, and a fixed volume to the capillary bed. In the brain, 

capillary density varies greatly among different tissue types. For the purpose of 

10 obtaining a robust model for all types pf brain tissue with the same choice of basic 

model parameters, the microvascular volume was kept as a free parameter. In terms 

of vascular transport, the models are similar^ except that in the model presented here, 

dispersion takes place only in the large* vessel, whereas small vessel dispersion fs 

I 

accounted for by the flow heterogeneity PDF. Despite the simplification of the model. 

»* 

15 It is likely that, by limiting the number. of vascular elements and thereby the number of 
free parameters, the stability of the m6del to noise has increased significantly, as 
outlined in the analysis above. The rob^ustness of the model to experimental noise Is 

Imperative to ultimately study CBF at high spatial resolution. Kroll at a/, assumed that 

i 

only small vessel tracer levels are observed by the residue detection (Kroll 1996). 
20 Here, this assumption is further Justified by the Inherent sensitivity of susceptibility 
contrast to microvessels (Fisel 1991; l!velsskoff 1994; Boxerman 1995). 

i . 

Vascular transport and dispersion terrn 

As discussed elsewhere t0stergaard 1d96b), nonparametric deconvolution approaches 

25 do not allow separation of macrovascCrlar transport end microvascular retention. 

Therefore* modelling vascular retention' Is necessary, whenever vascular transport 

significantly changes the input bolus s.^ape upstream of the AIF measurement site. 

The vascular transport operator modelfused in the present model is generally accepted 

for modelling normal major vessel transport. The advantage of including vascular 

30 transport In the kinetic modelling is de)monstrated in Rgure 6: The fact that tracer 

arrives earlier in tissue whh high flow fate due to a faster feeding vessel transit le 

accounted for by the model approach,! yjniike the modeNfree SVD approach. The 

difficulty In Introducing this operator l/es mainly in the fact that the transfer functions 

. of the major vessel and the microvascular network are very similar. This, In turn, leads 

I 
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to the difficulty in separating vascular v|lume and flow, as observed in the simulations 
(Rgure 9). The analysis presented here elso suggest that this effect contributes 
significantly to the uncertainty of flow &timates, even at nnodest noise levels^ just at 
it may play a role in the dependency upon initial conditions of fitted flow rates. These 
5 findings are In agreement with those prlviously reported by Kroll et al (Kroll 1996). In 
cerebrovascular diseases, blood may paps through stenoses with marked turbulence, 
or through irregular collateral paths upstream of the arterial sampling site. In these 
cases, the vascular operator may not bJ adequate, and the vascular transport furictlon 
should ideally be measured independenily. Indeed, novel MR techniques detecting the 



10 inflow of spin fabelfed arterial blood to 
information (Wong 1997), Such indepe 



given brain region may ultimately provide this 
dent measurements would serve to avoid the 



interdependence of flow and vascular v jJumes using vascular models, or alternatively 
allow application of the model-free SVC approach. In severe cases, however, vascular 
dispersion may dominate total tracer retention. In which case estimates of flow rate 
15 by residue detection of an intravasculai tracer become uncertain {(^stergaard 1996b; 
Kroll 1996). 

Flow heterogeneity 

Reports of cerebral flow heterogeneity |ire mainly based on invasive measurements of 
20 plasma or red blood cell velocities In ra ;s. Abounader ef a/, utilized bolus Injection of a - 
plasma marker followed by decapitatioi i, deriving plasma flow velocity from capillary 
filling in histology (Abounader 1995). Hudetz et determined the frequency of red 
blood cell IRBC) velocities in cepillarleslusing intravital video microscopy (Hudetz 
1997). Both found the distribution of b ood elements to be very heterogeneous, with a 



25 right-skewed shape. Figure 1 0 shows 
here. The data of Hudetz et oi (Hudetz 



heir findings along with the PDF determined 
1997) were normalised as the relative flow PDF 



(See Theory, Equation 9). From the dai a of Abounader et aL (Abounader 1995), a 
normo-capnlc data set was chosen, an J axes scaled to facilitate comparison with the 
other curves (plasma flow units did no| allow direct normalisation). It should be 
30 noticed that these measurements are i|ot directly comparable. Rrst, plasma and RBCs 
follow different paths through the capi lary network, so the measurements presented 
here should be more comparable to thi .plasma velocity measurements. Secondly, the 
flow distribution curve derived in this s tudy assumes equal capillary lengths. Although 
the relationship between capillary plasLa flow and flow velocity is a complex function 
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of capillary length and architecture, a finite distribution of capillary lengths Is likely to 
result in blood vetocities being more jlilspersed than relative capillary flows (As seen 
when comparing the PDF presented j[»erfe with the RBC velocity study). The similarity 
of the present measorements with tJiesp independent, invasive nnethods lends hope to 



5 the use of this approach in describlnj 



the heterogeneity of microvascular fiio 

found that plasma flow became morjji 

10 for red blood cell velocities was repfirti 

I 



srmal microvascular dynamics. In altered 



physiological states or disease, the f ow heterogeneity may not be as constant among 
tissue types as found In this study, if bounader et (Abounader 1995) determined 



at different degrees of hypercapnia, and 
mogenous at higher flows. A similar finding 
by Hudetz eta/. (Hudetz 1997). This could 
Ion should therefore be exercised jn choosing 



be the case in disease as well, and | au 
flow heterogeneity PDF for use with vafecular models In these cases. The model-free 
approach to determine flow hBterog|j n« ity PDF presented here may provide Insight 
with respect to the distribution of re; at ve flows in these cases. Using this approach, it 
15 should be kept in mind solving Equ J loi 1 belongs to a class of so-called Inverse 

problem, meaning that any noise in jrie \sured tissue concentrations may lead to large 
changes In the resulting residue fun^tic n. As suppression of noise Inevitably causes 
loss of underlying Information, the ^Vej deconvolutlon therefore may not yield the 
exact shape of the underlying re8id| s function. Likewise, one may not be' able to 
20 distinguish slightly different flow heji :>geneity PDF based on noisy measured 
concentration time curves. Although tl-je flow heterogeneity PDF and derived flow 
rates were found In agreement with'| ride pendent findings, it is therefore important to 
further validate the approach preserjleJ here. 



25 • Utility of vascular models 

Although-cerebral-blood-flow-itself-L 

heterogeneity of microvascular flb\Af3Brb transit times described here may be a more 



Important determinant of cerebral naJaTjboltsm. As discussed by Kuschinsky et a/,, the 



degree of heterogeneity among cjapnlaf / paths determines the net caplllary-to-tissue 
30 concentration gradients necessary tf ic -ive delivery of nutrients (Kuschiniky et al.. 



i-important-index-cf-braln-functior^thft- 



1992). Indeed, regulation 



Bcessaiy lb ii 
of capillajy ' 



f ow heterogeneity may play a major role in the 



brains ability to Increase e.g. oxygelrB d 3livery to meet cellular metabolic demands 

(Kuschinsky at al., 1992). This Issiil c m be addresses In great detail by combining 

• in ' ' 

flow heterogeneity measurements yjltll spatially distributed models of oxygen 
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exchange (U 1997). The analysis ^riented here, combined with these models, may 
ultimately lead to a more extcnslv;d |id6rstanding of, for example, the fundamental 
limitations of oxygen delivery in sjro|e. where the survival of tissue may partially 
depend on the ability to Increase i,x4en extraction by an Increased mean transit time. 
5 Also, modeinng the exact relatlon4h| between cellular oxygen consumption and 
vascular oxygen levels may faciirtktla quantitative metabolic interpretation of 
deoxyhemofllobin concentration cV,a^ges observed by functional magnetic resonance 
imaging (Kwong 1992). 



If 
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EXAMPLE 2. 

MAGNETIC RESONANCE (MAGING 
DEMONSTRATE HIGH RISK OF IN 



5 SUMMARY OF THE EXAMPLE 
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MEASUREMENTS OF FLOW HETEROGENEITY 

I 

F^RCTION IN ACUTE STROKE 



i 



The ability of brain tissue to surviyL ischemic episodes is believed to be related to 
changes in the dynamics of capillary blood flow. Using a novel magnetic resonance 
Imaging based method to measur* .mlcrovjascular transit time dynamics, the flow 
heterogeneity and tissue mean plaisma transit times was examined in 1 1 patients 

10 presenting with acute (<12 hour; rafter symptom onset) stroke. In normal brain tissue, 
the distribution of relative flows vr'as found to be markedly skewed towards high 
capillary flow velocities. Within regions of decreased cerebral blood flow, plasma 
mean transit times were prolonged. Furthermore, subregions were identified with 
significant loss of the high-flow c imponent of the flow distribution, thereby causing 

15 increased homogeneity of flow velocities.' These findings are In agreement with 
independent, invasive measurements of flow heterogeneity in states of decreased 
perfusion pressure In animal models. In p&rametric maps quantifying the acute 
deviation of flow heterogeneity frpm that' of normal tissue, areas of extreme 
homogenisatlon of capillary flow a predicted final infarct size on follow-up scans In 10 

20 out of 1 1 patients. Row heterogm;»eity and plasma mean transit time can be rapidly 
assessed as part of a routine clinfcal maghetlc resonance examination, and may 
provide a tool for individual planntaig of stroke treatment, as well as in targeting antf 
evaluating emerging therapeutic strategiels. 



25 BACKGROUND 



Acute stroke is the third leading cause of death, and the leading cause of adult 



disability. Emerging therapeutic strategies seek to minimise the progression of tissue 
damage in the acute phase of th<|^diseasi. Methods to rapidly assess the severity and 
later progression of acute stroke Ih indivi Jual patients are therefore highly desirable to 
30 plan individual treatment, as wellfas to evaluate novel therapeutic strategies. 

In acute cerebral Ischemia, deliva^ of nutrients is severely compromised, and tissue 
survival therefore depends on tissjje regi latory mechanisms to meet metabolic needs. 
Studies using positron emission t ojmogra jhy (PET) have shown that mean blood transit 
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time (MTT) and, with further dropjin cerebral perfusion pressure, also the oxygen 
extraction fraction (OEF) are mcr^sed inYtis&ue at risk' of infarction (Gibbs 1984; 
Baron 1981). Although the relaticliship bltween prolonged blood WTTT and OEF 
remains unclear, both phenomenJlare believed to reflect underlying regulatory 
mechanisms attempting to compepsate for a decrease in perfusion pressure. 



. One mechanism for vasoregulatorfy contrij^l Is believed to be the ability to alter the 
heterogene/ty of blood transit times and thereby the mean capillary concentration of 
substances diffusing from blood to tjssue[ (Kuschinsky and Paulson 1992). Animal 



decreajsed flow heterogeneity during whisker-barrel 
1 996) J indicating that this may ba the mechanism 
Ing ability to meet Increased metabolic needs during 



10 experiments in rats have revealet 
stImulation(Vogel and Kuschinsky 
underlying the normal brain's strf 
functional activation. FurthermonS, Hudeifz Bt ah demonstrated that graded decrease in 
perfusion pressure causes progressive lojs of high-flow components, thereby 

15 decreasing total flow heterogenenyCHudetz 1996). Decreasing flow heterogeneity 
therefore seems to play a cruclafjrole in maintaining sufficient concentration gradients 
to drive diffusion of nutrients sucp as ox rgen from blood into the cells (Kuschinsky 
and Paulson 1 992), This suggests that bj ood mean transit time and the degree of flow 
heterogeneity are important indic, ss to assess and further understand the ability of the 

20 brain to survive Ischemic episode \, 

The heterogeneity of flows in notmal volunteers by magnetic resonance Imaging (MRI) 
residue detection has recently be studied {0stergaard 1999). It was found that the 



probability density function |PDFj 
25 and among normal volunteers. In 



of relajive flows was remarkably constant within 
this study, magnetic resonance residue detection 



was used to study plasma mean!|ransii times as well as flow heterogeneity patterns in 



patlents presenting with acute st'oke 
with later neuronal death by corr Daring J; 
follow-up MRI or computed tome jraphy 



FCrthermore, these findings were correlated 
litial diffusion weighted imaging (DWI) with 
(CT). 



30 



The present findings show that f bw heterogeneity 
detectable by Invasive microscopy in an 
examination as part of routine MkI of acbte 
of heterogeneity change relative.^to nomjel 



changes, previously only 
Ttal models^ can be assessed by a 2-m1nute 
stroke patients. Furthermore, the degree 
tissue Is a powerful predictor of later 
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neuronal death, suggesting flow heterog^peity may provide an important diagnostic 
tooi in stroke patient management. 



MATERIALS AND METHODS 



Patient data 

All patients were treated wkh best mcdicjal management, but did not receive tPA or 
other thrombolytic iraeatment. DWI and CBF results for these patients are reported 



earIFer(Sorensen 1 998). 
10 Imaging was performed on a GE Signa 1:5 T imager {General Electric, Waukesha, Wl) 



retrofitted for EPI capabilities (Instascan, 
MRI Perfusion Protocol. Determination of 



Advanced NMR Systems, Wilmington, MA). 
CBV, CBF, MTT, and flow heteroQcnelty. 



Perfusion imaging was performed using ^In echo (SE) or gradient echo (GE), echo 
15 planar Imaging (EPI) with a time of repeti :ron ^R) of 1 .5 seconds, and a time of echo 
(TE) of 100 ms (50 ms for GE EPI). ThelsMce thickness was B mm with an In-plane 
resolution of 1 .56 mm by 1 .56 mm In a |40 by 20 cm field-of-vlew (FOV), A total of 
52 images were acquired, starting 15 seionds before i.v. Injection of 0.2 (SE-EPI) or 
0.1 (GE-EPI) mmol/kg Gd-based contrast jagent. Intravascular contrast agent 
20 concentrations, C(t), were quantified assuming a linear relationship between 

concentration and change in transverse ^jslaxatlon rate, AR^ (Villringer 1988; Weisskoff 



1994). The shape of the arterial input function (A)F) was determined from feeding 
anerlal branches, identified in the fmageiltce as pixels displaying early concentration 
Increase after contrast injection (Porkka ;i991 ). The tissue residue function (or Impulse 
25 response function) was calculated by deconvolving the tissue concentration time 



curve by the AIF, using singular value decomposition (SVD) (0stergaard 1996a; 
0stergaard 1996b). CBF was determlnec as the height of the deconvolved tissue 



curve. CBV was determined by the area 



jnder the tissue concentration time curve, as 



previously described (Rosen 1 990), and ,the plasma mean transit time (MTT) formed as 
30 the ratio CBV/CBF{Stewart 1 894). FInalljf, the distribution of tissue transit times //) 
each Imaging vox6/ was determined as t/na slope of the residue function, and by 
assuming equal lengths of capillary paths, the corresponding probability density 
function (PDR of relative flows waa determined by means of the central volume 
theorem (Stewart 1894; 0stergaard 1996). In order to quantify and compare the 
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deviatfon of the experiment ally determinld PDF from that found in normal brafn, a 

Kolmogorov-Smirnov test was performed! comparing the flow PDF in a given pixel to 

Iff 

that prevfously detemiined In normal tissue (See Results, Figures 11a, 11b and 11c) 
(Press 1992; 0stcrflaard 1999). The corr|spondtng p-value (Null-hypothesis: flow 
heterogeneity distribution is equal to thatjof normal tissue) was considered statistically 
significant if p<0.01 (Without Bonferonljborrection). 



Initial and final Infarct sizes. 
At the initial scan, infarct size was asse^ed by diffusion weighted imaging (DWI) 

10 (Moseley 1 990), acquired using EPI and j|-values ranging fronn 892-1 221 s/mm*. 
Measurements ware performed in 17-20jMces to cover the whole brain. Final Inferct 
size was assessed from DWI images acqiired 2-5 days after the Infarct, from or 
FLAIR MR images acquired at least 5 davjs after the infarct, or from CT images 
acquired more than 5 days after the infatet if no MRI was available. 

15 Volumetric Analysis || 

Using a semi-automatic Image analysis software package (ALICE, Hayden Image 
Processing Group, Boulder, CO) initial arjl final infarct volumes were calculated based 
on Initial DWI and follow up studies, restfectrvely. In addition, volumes of prolonged 



MTT and abnormal p (p<0.01) were de 
20 combined DWI and MTT, respectively p 



RESULTS 

Eleven hyperacute stroke patients, 7 ma 
25 BO) years, were examined within 12 hojl 
showed-diffusion-abnormalltles-on-the In 



jned, In order to assess the ability of 
laps to predict final infarct volumes. 



and 4 female, with mean age 61 (range 33- 
s of symptom onset (Table 3). All patients 
tial DWI {MoseleyH 990)i-consrstent-wIth— 



Ischemic neuronal death prior to. the Inrtii tl scan. Furthermore, all patients showed 
greater volumes of decreased carebrel b'i jod flow (CBFJ and/or increased MTT In the 
hemisphere of infarction. • 



30 
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Patient 


Age/sex 


Hours 


1 


33 M 


4.5 


2 


7^F 


3 


3 


65 M 


2 


4 


53 M 


3 


5 


72 F 


2.5 


6 


64 F 


5.5 


7 


79 M 


7,0 


8 


BO F 


1V0 


9 


45 M 
65 M 


4.0 


10 


11.5 


11 


45 M 


6.5 



15 



Age, sex, tim© of initial MRI tor j1 1 acute 

ii 

spontaneous reperfusion of the occluded 



1042 



stroke patients. Patient 9 showed 
vessel on the follow-up MR angiography. 



20 Flow heterogeneity findings 

Outside volumes of profongad MTT, the ^hepe of the tissue flow PDF was similar to 
that previously found In normal k^olunteer namely a right-skewed distribution with a 
distinct distribution of high flov\j! rates. Im^lde volumes of prolonged mean transit 
tinnes, the shape of the flow PDF was eif^er tike that of nomnal tissue, or showed 
25 distinctive loss of the high-flowiportion of the PDF; To Illustrate the first type. Figures 
11a, lib, and 11c show a typical pattern in patient 6. This patient, a 64 year old 



3ici 

i 



female, was exannlned 5.5 houTjs after or set of left leg weakness, and showed 
prolonged MTT corresponding to the ant< rior cerebral artery territory (Figure 11a). 



plots for normal brain tissue as well as two 



Figure 1 lb shows the flow heterogenefty 
30 regions of prolonged mean tran iit time [A seas are Indicated on the MTT map with 
numbers corresponding to the F OF curve s). The flow PDF In normal tissue was 
markedly right-skewed, and matched the 
(0stergaard 1999). The volumes of incrc 



shape previously found In normal volunteers 
ased MTT displayed PDFs with a more 



symmetric shape, with a tender cy to loobe the high-flow population found in normal 
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tissue. The degree of synnmetry varidii within the volume of increased MTT. The 
deviation from the normal PDF wlas sTubsequently quantified by a Kolmogorov-Smlrnov 
test, yielding the probability p that the curve belong to the distribution of relative 
flows of normal tissue (from ref i0stlrgaard 1 999)). In figure 11c, areas with large 
deviations of the PDF from that of normal tissue {p<0.01) arc shown by a colour 
coded overlay of p onto the acutje'CaF map. Based on the previous experience 
(0stergaard 1999) a significance fevel of p<0.05 displays few PDF abnormalities in 
normal tissue, except in major vessels. Therefore, p<0,01 was chosen to highlight 
highly significant deviations froni normal flow heterogeneity PDF. 

MTT, flow heterogeneity and later inifarction 
In 8 out of 1 1 patients, comparison bf initial DWI Images with the follow-up study 
showed that lesion size had increased in between the Initial and fot1ow*up scans. In ail 
8 cases, neuronal death had occWrep within the region initially displaying increased 
15 MTT, Figures 12a and 12b show thiS correlation In patient 11, a 45-ycar-oId male 6.5 
hours after onset of symptoms. On the Initial DWI (first row), cell death is localised to 
deep gray matter, whereas the acute MTT maps (second row) show prolonged MTT 
corresponding to the whole middle carebral artery territory. The p-map shows highly 
significant deviations from the normal flow PDF in anterior and posterior sub-reglons. 
20 These subregions corresponded , well with tissue that later infarcted, displayed as 

bright regions In the 2-month follow up FLAIR MRI images (bonom row). In Rgure 13, 
final infarct volumes are compared tjo the initial abnormalities of DWI + MTT and 
DWI+p maps, respectively. Notice we striking ability of combined initial DWI and p- 
maps to predict final Infarct size'. ;In m 0 out of 1 1 patients, p maps (excluding vessels 
25 and volumes with preserved, high flow component, see below) corresponded well 
^wlthjtheJInaiinfarc.tJBgur_esJij4:€(_andJj4^ 



Notice that the initial lesion (top rovCj and extent of MTT prolongation (second row) 
are similar to that observed In Jrigjurls 12a and 12b. Although the CBF was 
significantly decreased (Third r<^, gray scale Image), there were only small 



30 abnormalities In the flow heterolgeneitles (indicated by coloured areas on overlay on 

• ' u 

the CBF map). The small spots correspond to vessels, having a high degree of flow 

vQS similar to that observed on the initial DWI, 
again have served as a predictor of final 



heterogeneity. The final Infarct sizef 
indicating that the heterogeneity ma 
autcomft. ^^ 



i : 
I : 
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To demonstrate the predictive value of the p-maps In cases of white matter Ischemia, 
Rgures 1 5a and 1 5b show mapsifrom patient 8, a 80 year old female. Small Infarcted 
areas are seen on the initial DWHrnages {Top row). The MTT Image shows 
5 abnormalities extending into whl^e matter around the Initial lesions. Again^ flow 
heterogeneity changes are observed in smaller subregions, corresponding well to the 
regions that went on to infarction. Wotice symmetrically located areas with low p 
values, corresponding to vessels! ! 



10 Artifacts in p-maps 

In one case, (patient 8), p maps Linderestlmated the final infarct siM. Figure 16 shows 
one slice from this patient, displaying areas with p<0.1. Notice high-intensity areas 
correspond to areas that later Intarcted, whereas a number of areas showed unspecific 
p increase. Separate evaluation of areas with p<0.01 would In this patient causa 

15 underestimation of final Infarct sSze.' Interestingly, follow-up MR angiography In this 
patient demonstrated spontaneojjs reperfusion between initial and follow-up scans. In 
the p-maps, small areas of low f|-values were In some cases observed at the location 
of major vessels (See Figure 12b and 15b), due to the homogenous flow pattern in 
vessels relative to that In tissue.! These areas were not included when defining areas 

20 of abnormal tissue flow heterog(tneity for comparison with follow-up studies. In 
patients 3 and 4, areas unrelated to'major vessels showed low p-values in a single 
slice, whereas adjacent tissue Inf neighbouring slices showed no abnormalities. 
Analysts of the flow heterogenerty PDF in these single slices revealed a high-flow 
distribution similar to that of normal tissue, whereas the low-flow component showed 

25 flow components down to zero {unlike the relatively sharp cut-off at 0.5 observed in 



normal tissue - See Figure lib) 
AIF relative to the tissue. Based 
PDF observed in adjacent tissue 
when comparing p-maps with 
30 discussed. 



DISCUSSION 



This Is Interpreted as being due to dispersion of the 
on the preserved high-flow component and the normal 
In neighbouring slices, these areas were not Included 
follovtf-up images. Below these phenomena are further 



The study confirms the report by Hadetz ef a/in animals that decreased perfusion 
pressure {CBF:CBV ratio) is associated with progressive loss of high-flow 
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components(Hudetz 1996). The ^udy extends these findings by documenting loss of 
flow heterogeneity In human acute stroke, and good agreement between 
heterogeneity changes in early Isdhemia and eventual tissue Infarction, The findings 
presented hence strongly support the hypothesis by HudeU et a/ that gradual loss of 
5 the high flow component of the f ow heterogeneity PDF heralds local loss of functional 
reserve capacity, and thereby neuronal death (Hudetz 1996). The finding that p-maps 

certainty iniuntreated patients suggest that MR 
heterogeneity measurements ma> prove useful for individual planning of patient 
management, as well as for evalijatlon of new therapeutic approaches in smaller 
10 patient populations. 



The hypothesis of heterogeneity changes being the driving force in regulating oxygen 
delivery to tissue (Kuschinsky and Paulsoh 1992; Vogel and Kuschinsky 1996) 
suggest a close relationship betxAjeen the findings presented here and the OEF increase 

15 observed by PET in tissue at high risk of subsequent infarction{Powcrs 1991; Wise 
1983). Indeed, qualitative analysljs of the kinetics of oxygen delivery show that one 
should expect reduced heterogerlelty of blood flow to produce an increased flow of 
oxygen into the tissue In states c f decreased. flow, as Illustrated In Rgure 17. The 
curve Is a plot of the oxygen flow' Into tissue versus blood flow. From the convex 

20 shape of the curve h can be seen .that oxygen flow into the tissue at a given mean 

blood flow Is greater when blood flow is homogenous than when blood flow is more 

heterogeneous: Oxygen flow Into tissue versus blood flow F^. related through the 

I ! 

the permeability to oxygen times the surface area of 
983).- If all blood flow is at the normal mean flow, 
given by the height of A. If part of the flow is at 
jhtings fo rnaintain the semeTTiBan^lowrth^o^cygen — 



equation F,.(l-e ''^j), where PS is 
capillaries (Renkln 1959; Crone 

25 

f norm* 

the oxygen flow into tissue 
and-the rest-at-fni^wlth-wei 



and fto« changed to maintain the 



flow Into the tissue will be the haFght of D. Notice, as mean flow is reduced and f^ 



same CMROj, both f^,^ and approach f^^, thereby 



decreasing the degree of heteroji^nelty. Notice decreased flow heterogeneity with 
30 constant flow wHl increase oxygfen delivery, in parallel with the findings of Vogel ef a/. 
who found decreased heterogeneity In states of functional acUvatlon and thereby 
increased oxygen metabolifim{VAgel and Ku^hinsky 1996). 
The observed shifts toward a hcmogenous flow distribution may therefore signal 
Increased utilization of metabolic regulatory capacity, explaining the risk of neuronal 
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death observed in these regions c 
method of choice to demonstrate 
disease. However, future studies 
heteroQeneJty meesurements end 
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^extreme flow homogenlsation. PET is today the 
rinetaboilc reserve capacity in cerebrovascular 
sjhould focus on the relationship between MR flow 
OEF measured by PET to further explore this 



5 couplinfl of microvascular dynamicjs and metabolism. 

! 

Neuronat death was localised in ereas Initially displaying prolonged plasma mean- 
transit times. This is agreement with the previous experience using this technique 
(Sorensen 1998), as well as studies using PET (Heiss 1994; Baron 1981) and single 
10 photon emission computed tomography (SPECT) {Buell 1 988). Although the CBVrCBF 
ratio {i.e. 1/MTT) depends linearilyjon the cerebral perfusion pressure over a range of 
values(Schumann 1998), this dependence Is likely to be lost when maximum 



vasodilation is reached at low pres3ures(Powers 1991). The MTT prolongation may 



therefore not be directly related to 



15 hence risk of infarction. The find 
early sign of decreased perfusior!) 
may still suffice to ensure tissues 



Atthough the high flow compone 
20 component of the flow heterogei 



the severity of the perfusion pressure drop and 



rigs support, however, that prolonged MTT is an 

>ressure, at a stage where regulatory mechanisms 
stirvlval. 



njc seem crucial to tissue survival, the low flow 
ity PDF may also prove useful In planning 



therapeutic approaches. Intra vitajljmlcroscopy studies suggest that maintaining 
capillary flow velocities above a.'tixed, lower limit Is essential to avoid white blood cell 
plugging of capillaries fHudetz 1^6; Yamakawa 1987). The distributions of absolute 
flows in single pixels may prove | useful in assessing of leukocyte adhesion prior to 
25 therapeutic attempts to reperfus ^| tissue. 



In patients with cerebrovascularj disease, the AlF may undergo dispersion and delays 
upstream of site of measuremen t| possibly causing overestimation of MTT Idstergaard 
1996a; Ostergaard I9g6b). Thi^ bias was reduced by choosing AlFs In the vascular 
30 territory affected by the vascula jbcclusi6n. Furthermore, dispersion of the AlF will 
tend to broaden the flow PDF, nWefore, the effects of dispersion counteract the 
observed homogenlsation of floy 'felemerfts. In determining the flow PDF, high p 
values were observed near ve5^^l5 (and therefore easily identifiable on the 

:i^jor vessel flow Inherently homogenous. Probability 



accompanying CBV maps), as rn 
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maps should therefore be careful 
signs of vessel dispersion in the .^1 



inherent microvascular weighting of these 
Weisskoff 1 994). 



Given these precautions^ the fin 
based assessment of flow heter 
metabolic reserve capacity in per^i 
10 angiography, DWI, determlnatiorj 
can be performed In roughiy 20 irn 



SI047 



y^inspecced for vessels on CBV maps^ as well as 
F shape In a given region, SE-EPl Images are 



particularly suited for this type o 'analysis, as large vessels are suppressed due to the 



images (Fisel 1991; Boxerman 1995; 



ngs presented here Indicate that magnetic resonance 
sneity provides a powerful too! to study residual 
nf arct tissue. Combined conventional MRI, MR 
])f flow heterogeneity and plasma mean transit timo 
jnutes on most clinical MR systems. Unlike PET and 
SPECT, examinations can theref{i(e be pJrformed within the short time-window where 

er to f revent the progression of neuronal death, 
hete ogeneity may In the future serve to guide 
1 nb pointjto tissue that may serve as target for novel 



treatment should be initiated In ^ 
.Presence of tissue with loss of 
15 Individual patient management, 
therapeutic approaches. 
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CLAIMS 



1 , Method for determining haemc}dyn|tnic ndices of an organ or of a part of tissue of 
a mammal including 



a) determining a time series of tc^og 
tissue during and after a bolus inject 
being substantially Intravascular 



10 b) determining a time series of 

concentration of the tracer in art4'rles|bf xh 
tomographic data. 



cc^nce^tration data being indicative of the 

e organ or tissue from the time series of 



c) determining a residue function 
15 of the time series of tomograph 

d) determining a distribution of trjanslj 

2. Method according to claim 1 /|s)vhe! 
20 normalised haemodynamic index is d 

the index being normalised by th 

3. Method according to clafm 2, 
a quantitative haemodynamic pa 



25 



4. Method according to claim 3; 



12)053 



data pertaining to the organ or part of 
a tracer dose to said mammal, the tracer 



organ or of the part of tissue by deconvolution 
:wltn the time series of concentration data, and 

^imes from the slope of the residue function. 



ein 



probability density function (PDR of a 
iermtned from the distribution of transit times, 
the integral of said Index. 



vafije of 

1 



I 

whe'irein 

a git 

■pm^er 



Qt least one of the haemodynamic indices Is 
obtained from the PDF. 



li 



vh^'ein IX least one of the at least one parameter Is 



obtain from comparison of the dftcnfiined 
PDF, 



PDF and a previously determined reference 



30 5. Method according to claim 4/ 
Kolmogorov Smimov test. 



Whd^^ein tSe parameter Is obtained by use of the 



Hi 



6. Method according to claim 3 ^nd Ipmp 



sing the steps of 
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determining the Impulse response fur«:tion of the organ or of the part of tissue by 
deconvolution of the time series pi tomographic data with the time series of 
concentration data, 

5 

determining the relative tissue fltJw frpm the Impulse response function of the organ or 
of the part of tissue, 

normalising said time series of concen tration data with the integral of said time series 
10 of concentration data with respept td time. 



volume, of the organ or part of tissue 
series of normalised concentration da 



15 



determining the normalised relative tiiisue flow, respectively the normalised blood 



by use of the relative tissue flow and the time 
a, and 



converting said normalised relatii-e tissue flow, respectively normalised blood volume, 
to an absolute value for the tissue flol^v, respectively the blood volume, by means of a 
previously determined conversion facjior. 



1 



1 



blood volume. 



I 



30 



eter being of metabolic significance and 



20 the quantitative haemodynamic parar 

determined from the PDF and the absolute tissue flow (FJ, respectively the absolute 



7. Method according to claim 3 and claim 6, wherein a parameter (E) significant for 

25 the local extraction of a substance Isjdetermined, the method further comprising the 

following steps: j 

calculating the relative fldw heierogGnelty (w(f)1 as a function of^he relatlvo~ 

flow (f) from the distribution of transil 

estimating a value (P) for the bcal;capil[ary permeability, 

estimating a value (S) for the iDcal; capillary surface area, 

calculating said paramete (E);Ls the integral value of the relative flow 

heterogeneity {w(f)) multiplied by one minus the natural exponential function of the 

I 

negative ratio between ) 
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i) the product of the locel capillary permeability (P) and the local cepllfary surface 
area (S), and 

ii) the product of the relative flow (f) and the absolute tissue flow (FJ 
with respect to the relative flow (f). 

6 

8. Method according to claim 6 or 7, wherein the normalised relative tissue flow, 
respectively the blood volume. Is also normalised with the ratio between body weight 
of the individual mammal and the injected tracer dose. 

10 9. Method according to any of claims 6-8, wherein the previously determined 
conversion factor is in general applicable for the present method to members of a 
mammalian specie. 

10. Method according to any of claims 6-9, wherein the previously determined 

15 conversion factor is In general applrcable for the present method to an organ or tissue 
of the mammalian specie. 

1 1 . Method according to any of claims 6-1 0, wherein the previously determined 
conversion factor Is a constant factor applicable for the present method for any organ 

20 or any part of tissue of the mammalian specie. 

12. Method according to any of claims 6-1 ^, wherein the previously determined 
conversion factor is a constant factor applicable for all of cerebral tissue of the 
mammalian specie. 

25 

13. Method according to any of the preceding claims, wherein the tomographic data 
are obtained by means of magnetic resonance imaging. 

14. Method according to any of the preceding claims, wherein the tissue is cerebral 
30 tissue. ' 

15. Method according to claim 14, wherein the tomographic date are obtained by 
means of susceptibility contrast magnetic resonance imaging. 
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I 

I 

1 6. Method according to any of the preceding claims, wherein the tomographic data 
comprfse Information pertaining to subreglons of sections of the organ or part of 
tissue and the heemodynamic Indices are determined for at least a substantial part of 
said subregions. 

17, Method according to claim 3 and 16, wherein quantitdti\/e haemodynamtc 
parameters are represented as Images subdivided into a plurality of pixels each 
representing a quantitative haerriodynamlc parameter pertaining to ona of said 
subregions. 



10 



15 



18. A system for processing of time series of tomographic data pertaining to an organ 
or a part of tissue according to any of claims 1-17, said system residing on a 
computer having means for producing an output representative of at least some of the 
determined haemodynamic Indices. 

19. Method for evaluating the efficacy of a drug or a substance on an organ or on a 
part of tissue of a mammal by rneans ot haemodynamic indices of said organ or of 
said part of tissue obtained by a method according to any of claims 1*17. 

20 20. Method according to claim 19, wharein a system according to claim 18 Is used. 

21. Method for obtaining Information of the likelihood of recovery of an organ or part 
of tissue in a living mammal upon or during a period of Insufficient vascular supply of 
said organ or of said part of tissue in the mammal comprising determining 

25 haemodynamic Indices accordinis to clalma 1-17. 

22. Method for obtaining information relevant for discrimination between relevant 
therapy of an organ or pan of tissue In a living mammal upon or a period of 
Insufficient vascular supply of said organ or of said part of tissue In the mammal 

30 comprising determining haemodynamic' indices according to claims 1-17. 

j 

23. Use of Information obtalne<| by usef of the method according to any of claims 1-17 
or any of claims 19-22 for preparing a reference table for use In discrimination of a 
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treatment schedule for an individual mammal or group of mammals for which 
information have been obtained In a manner similar to said information. 
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